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Abstract: We consider the maximal cut of a three- loop four point function with massless 
kinematics. By applying Grobner bases and primary decomposition we develop a method 
which extracts all ten propagator master integral coefficients for an arbitrary triple-box 
configuration via generalized unitarity cuts. As an example we present analytic results for 
the three loop triple-box contribution to gluon-gluon scattering in Yang-Mills with adjoint 
fermions and scalars in terms of three master integrals. 
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1 Introduction 

Precise predictions for cross sections in collider experiments require knowledge of scattering 
amplitudes to a high loop order. Computations of this type, particularly in QCD, are 
still unfortunately out of reach using known techniques. Unitarity cuts have shown to 
be an enormously powerful tool for super-symmetric gauge and gravity theories with full 
amplitudes now computed up to four loops [1, 2]. At three-loops, steps beyond maximally 
super-symmetric amplitudes have also been considered in the recent computation of the 
UV behaviour of graviton-graviton scattering in A/" = 4 supergravity [3]. For a recent 
review of state of the techniques see [4] and references therein. 

The success of unitarity [5, 6] and generalised unitarity [7] in automating multi-loop 
one-loop corrections has sparked some recent interest in exploring the application of inte- 
grand reduction (OPP-like [8]) methods at two- loops [9-14]. Very recently the possibility 
of using computational algebraic geometry to overcome the traditional bottlenecks in am- 
plitude computations has started to be explored [15, 16]. 

The traditional approach to a unitarity computation of a loop amplitude relies on the 
knowledge of a basis of known integral functions. At two loops it has been possible to 
derive sets of integration-by-parts identities [17, 18] that reduce amplitudes to a unitarity 
compatible basis [19, 20]. Using an integrand parametrisation constrained by the ideal 
generated from the propagators with help from Gram matrix identities we can also reduce 
a unitarity compatible form of the amplitude to master integrals using the well known 
Laporta algorithm [21] which is implemented in a number of public codes [22-25]. 

In this paper we consider the extension of the new integrand reduction techniques to 
three-loop amplitudes. We derive a complete reduction to ten-propagator master integrals 
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(MI) for the maximal cut of the three loop planar triple box with four massless external 
legs. As a first application of the technique we compute the ten-propagator MI coefficients 
for gluon scattering in Yang-Mills theory with adjoint fermions and scalars. The result 
applies in both M = 4,2,1 and super-symmetric theories. 

2 Reduction of the Massless Triple Box 

In this paper we will study the three- loop planar triple box topology, shown in figure 1, 
defined by: 



where the ten propagators {/j} are given by: 



(2.1) 



h = ki, 

h = h+Pi +P2, 
Ig = ki + ks. 



h = ki -pi, 
h = h -P4, 
ho = k2- k^. 



h 
h 



h -pi -P2, 



k = h+Pi +P2, 
h = h, 



(2.2) 



The external momenta, {pi}, and internal momenta, {h}, are considered massless. 
Though the integral needs to be dimensionally regulated in D = 4 — 2e dimensions we will 
only consider generalized unitarity cuts in four dimensions and therefore reconstruct only 
the leading term of the numerator function, Nni^, ki,Pj) = N{ki,pj) + 0{e). 

We also consider the external momenta to be outgoing with the Mandelstam variables 
defined by: 



(pi +P2)^ 



t = iPi+PAf 



-s-t. 



(2.3) 
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Figure 1. The momentum flow and propagator definitions for the three- loop planar triple box. 



- 2 - 



Following a similar approach taken at two loops in our previous work [11], we will 
proceed in three steps: firstly the on-shell constraints for the deca-cut will be solved. 
Secondly, we invert a linear system to map the polynomial of the deca-cut defined by the 
on-shell constraints onto the general integrand basis. Finally the integrand is reduced to 
master integrals (Mis) using additional integration-by-parts(IBP) relations. 

2.1 Solving The On-shell Constraints 

The parametrisation of the loop-momenta using two component Weyl spinors has been 
chosen as follows: 



h = xipi + X2P2 + X3 



-k = yiP3 + 2/2^4 + ya 



-k = Z1P2 + Z2P3 + Z3 



(23) {pi\r\P2] 

(13) 2 
(41) (P3|7^b4] 
(31) 2 
(34) (P2|7^|P3] 

(24) 2 



+ X4 



+ 2/4 



+ Zi 



(13) {P2\r\pi 

(23) 2 
(31) (P4|7^N 
(41) 2 

(24) {P3\1^\P2] 

(34) 2 



(2.4) 



The ten cut-constraints if = form a set of ideals. Following the recent method proposed by 
Zhang [15], this system can be reduced using primary decomposition to find 14 independent 
solutions. These solutions come in complex conjugate^ pairs and we label them 1 — 7 and 
1' — 7'. We find that all 14 solutions have the same dimension and therefore can be 
parametrised with two variables which we call ti and T2. For the solutions 1 — 7, explicit 
forms for the coefficients in (2.4) can be written as in tables 1 and 2. 
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Table 1. Values of the coefficients in I2 and —Iq in (2.4). 

It is worth mentioning that the parametrisation is Laurent-polynomial, so that no 
terms of the form appear. We will therefore be able to fit the integrand using an 
efficient discrete Fourier projection just as used succesfully at one- loop [26, 27]. 

The conjugate solutions can be easily constructed from the above expressions using. 
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(2.5) 



^Technically explicit complex conjugation is only valid for real external momenta. However, the solutions 
1' — 7' are also valid for complex external momenta. 
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Table 2. Values of the coefficients in -U in (2.4). Z4i = ^(l + ^jn + :^ 



Note that we have suppressed the functional dependence of the coefficients, e.g. x| = 
x|(Ti,r2). 

2.2 Fitting the Integrand Basis 

The space of loop momenta is spanned by three external momenta, which we choose to be 
{pi-iP2-,Pa\-, and an additional orthogonal direction oj which is defined by, 

^ = ^((2|3|l](l|7n2]-(l|3|2](2j7ni]). (2.6) 

The integrand numerator N , can be parametrised in terms of seven irreducible scalar 
products (ISPs) as follows: 

N = Y,(^a,...aAkl ■ PaT' ■ Pir ih ■ PATHh ■ PlT' 

x{h-ior'{k2-Lor^{h-ior\ (2.7) 

We can use renormalizability conditions to restrict the maximum powers of Oj and 
complete the reduction by using polynomial division with respect to a Grobner basis con- 
structed from the propagators constraints. This procedure is implemented in the Mathe- 
matica package BasisDet and we refer to [15] for full details of the method. The resultant 
integrand naturally splits into 199 spurious (S) and 199 non-spurious (NS) terms of the 
form: 

Tv^s = ■P4r{k2-pinh-PAr{h-pir^ 

X {cZ..a,0 + C^taAkl ' ^)(^2 ' ^)) , (2.8) 

= ^(fci ■ PAr{k2 ■ piTHh ■ p^nh ■ pir^ 

X {cl,...aAkl ■ ^) + cl,...aAk2 ' ^) + cl^.^aA^S ' ^)) ■ (2.9) 

The spurious terms will vanish after integration since they are linear in fcj • cj. 
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Inserting the expressions for the loop-momenta given in (2.4) into the integrand N = 
N^^ + , defines 14 Laurent-polynomials of ti and T2, 



N\s = J2dsijTl4, (2.10) 

ij 

where s denotes one of the on-shell solutions either 1 — 7 or 1' — 7'. In total we find 622 terms 
in these expansions which we collect into a vector, d. After collecting the 398 coefficients 
of (2.7) into a vector c we can define a Matrix M such that, 

d = Mc. (2.11) 

After inverting this matrix we can re-construct an arbitrary integrand by using the Lau- 
rent expansion of the product of eight tree-level amplitudes to extract the values of the 
coefficients, dsij. 

Though in principle we could invert the matrix in its full 622 x 398 form this task is 
quite complicated in practice. However, it is straightforward to separate the problem into 
two pieces corresponding the spurious and non-spurious terms. We can show that, 

k'{Tl,T2) ■p = k''{Ti,T2) -p, k'{Ti,T2) -UJ = -k''{Ti,T2) -UJ. (2.12) 

Therefore when we combine the paired solutions we find, 

A^l, + N\s' = N^^ls + N^^ls' = 2N^^\s, (2.13) 
A^l, - N\,, = N^\s - N^\s' = 2N^\,. (2.14) 

After defining two new vectors, d±, where 

4 = , (2.15) 

we find smaller matrices 311 x 199, M±, satisfying 

d+ = M+c^^ d_=M_c^. (2.16) 

In this form it was possible to invert the matrices analytically using standard computer 
algebra packages. 

2.3 Alternative Branch-by-branch Fitting 

In this section we describe a technique that allows each branch of the on-shell solutions to 
be considered as a separate linear system, a natural extension of the strategy taken at the 
end of the previous section when we separated spurious and non-spurious terms. This will 
be achieved by first partially fitting the integrand on the 14 solutions to get 14 polynomials. 
Then the 14 polynomials are combined together to get the full integrand, by a Grobner 
basis method. Some of the more mathematical details are described in the Appendix A. 

On each of the 14 solutions the integrand numerator A'^ can be reduced further to 
a polynomial Ni with much fewer monomials, by the multivariate division towards the 
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number oi terms 
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25 



Table 3. Number of terms for reduced integrand in each solution. 

Grobner basis of the corresponding branch. We note that Ni = A^j({ISPs}) is distinct from 
N\i = A^|j(ri,r2). The number of monomials of each Ni is listed in table 3. 

For example, divided by the Grobner basis of the second branch, the integrand is 
reduced to a linear combination of 59 terms, 

( A A 43 '^4 42 42 33 33 24 24 4 4 32 

y ,x y ,x y ,x z ,x y ,x z ,x y ,x z ,x y ,x z,x y,x z , 

322323 4 443 3 2222 3 3443 

X y ,x z ,x y ,xz ,xy ,x ,x z,x y,x z ,x y ,xz ,xy ,z ,y ,x , 

22 22332 22 

X z,x y,xz ,xy ,z ,y ,x ,xz,xy,z ,y ,x,z,y,l}, (2.17) 

where x = ki ■ p4, y = k2 ■ pi, z = ■ p^. The coefficients are determined by polynomial 
fitting techniques, at the second solution. 

d(2) = M(2)c(2), (2.18) 

where c^^^ is the list of the 59 coefficients and d^^^ contains the Laurent expansion coeffi- 
cients at the second solution. M^^^ is a 62 x 59 matrix, which is much smaller than M, so 
it is easy to compute c^^^ and the reduced integrand 

Once all 14 A^j's are obtained, we "merge" them together to get the integrand N. The 
step is equivalent to solving congruence equations in a polynomial ring. It is automatically 
done by a Macaulay2 program [28] , based on a computation using Grobner basis and ideal 
intersection. A more detailed mathematical description of the procedure is outlined in 
Appendix A. 

3 Results for gluon-gluon scattering in Yang Mills theories 

The starting expression for the computation is the product of tree-level amplitudes: 

= EE^/i-/.(^/'^^) 4°\-^2l^P2^^3:^) 

fi hi 

V /l('^V_/~'*3 _7~'^4 1^9 ^ /l{0)/'//l4 _7~h5 rhlO \ /l(0)/r/l5 A3 _]-hfj\ 

^ ^3 I '3,/i' %/2'*9,/4^ ^3 l'4,/2' '5,/3 '^lOJs^ ^3 l*5,/3'?^3 ' %h > 

V i'^'^V/'^e _]-h7\ /tCOV/ZiT _l-hs _]-hio\ A(^)(lh8 ihi _]-hg\ /-q 1 \ 
^^3 l''6,/3'-P4 ' '7,/3''^3 V''7,/3' '8,/2' 'lOJs / ^3 v'8,/2' 9,/4''' ^ ' 

where li and pj are defined in (2.2). {hi} are the internal helicity states and {Aj} are the 
helicities of the external gluons. There are 34 distinct configurations of the internal flavours 
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{fi} each associated with a number of fermion flavours, Uf and complex scalar flavours, n^. 
The explicit expressions for the flavour coefficients Ff^ j^{nf,ns) are given in Appendix 
B. 

Complete expressions for all integrand coefficients have been found as functions of 
and Us- Though each term is relatively compact, the full set of 398 coefficients makes 
the expression rather lengthy so we only present the result after further reduction to Mis. 
The full expressions for the integrands are available from the authors on request. We note 
however that all coefficients in the Laurent expansion (2.10) with \i\ or \j\ > 4 vanish for all 
Uf and Ug demonstrating the general renormalization conditions for this case are simpler 
than assumed when constructing the integrand in section 2.2. 

IBP relations generated with Reduze2 [24, 25] reduced the 199 non-vanishing integrals 
in the triple box integrand onto 3 Master Integrals, /io[l], Iio[{ki +^4)^] and /io[(A:3 — ^4)^]. 
The analytic expression for the scalar integral Iio [1] has been known already for some time 
[29]. We should be clear that to obtain the reduction of a full amplitude to all master 
integrals the procedure must be carried out in two steps. Firstly the full integrand, without 
applying IBPs, must be kept as a subtraction term for each lower propagator topology. 
Secondly IBP relations, again including all lower propagator Mis, should be applied to the 
complete amplitude. 

In order to map the basis of 7 ISPs of (2.8) onto the 15 propagator^ integral topology 
obeying IBP relations we must re-write the non-spurious term {ki-uj) {k2-uj) in terms of the 
reducible quantity /ci • /c2 . This is achieved straightforwardly using the Gram determinant 
identity generated by: 

de.^^'l'^f'Uo, (3.2) 

\12 4:ujk2j 

which is re- written as: 

{h ■ u;){k2 .u) = --+ '-{{h ■ 4) + {k2 ■ 1)) + -{k^ • A^s) + ^^(^i • A){k2 • 1). (3.3) 

4t ^ S S 

We note that this does not change the number of non-spurious integrals in the integrand 
basis. 

We write the 4-point 3-loop primitive amplitude for this ladder topology as: 

y4f^(1^^2^^3^^4^*) = Ci{l^\2^\3^^,4^^)Iw[l] 

+ C2{l^\2^\3^',4''')ho[{ki+PAf] 

+ C3 (1^^ , 2^2 ^ 3A3 ^ 4A4 ) P3 _ p4)2] + . . . (3.4) 

where the ellipses cover terms with less than ten propagators. We then define dimensionless 
coefficients by, 

Ck{l^\2^\3^^,4^^) = shA^^\l^\2^\3^^ ,4^^)C^^^^^'^*{s,t), (3.5) 

^For the purposes of the IBPs we consider both positive and negative powers of the propagators in the 
topology hence the 15 propagators includes both the 10 denominators of eq. (2.1) and the 5 non-spurious 
ISPs. 
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where A^^^ are the tree-level helicity amplitudes (expressions for the tree-level amplitudes 
are collected in Appendix B for convenience). After following the procedure of computing 
the Laurent series of (3.1), reconstructing the integrand and applying the IBPs we obtain 

the following results for the three independent helicity configurations. The turns 

out to be trivial, 

(7f-++(s,t) = -l, (3.6) 
(72--++(s,t) = 0, (3.7) 
C^-++{s,t) = Q. (3.8) 

The — h H — configuration is, 

C'f++-(s,t) = 

- 1 - (4 - n^)(3 - n,) + (4 - rif) 

-{l + ris- rif)^ {2t^ + list + lOs^) , (3.9) 

C-++-[s,t) = j{l-C^++-), (3.10) 
C'3~++-(s,t) = 0, (3.11) 
and finally the alternating — | h, 

(7r+-+(s,t) = 



St .sH"^ 



I + {A - nf)— - 2{l + Us - nj)- 



u 



, , , , sH{2t - s) 

+ (2(1 - 2n,) + ny) (4 - n^) > 

- (n,(3 - - 2(4 - n ff) ''^'" ' ^ (3.12) 



C2-+-+(s,t) 



s'^{2t - 



(2(1 -2n,)+n/)(4-n/)- 



^i■ 



4 



+ (n,(3 - usf - 2(4 - Uff Y^'' + (3.13) 

3s2(2t - s) 



+ {2{1 -2ns) +n f) {A- UfY- 



2^4 

/ / ^■? 3s(t^ - 4st S^) 

- (n;(3 - ?i,)2 - 2(4 - n^f)^ ^. (3.14) 

It is easy to check that these coefficients correctly reproduce the known result in A/" = 4 
super-symmetric Yang-Mills [30]. An expression valid for > super-symmetric gener- 
ators can be found by setting nj = M and = M — \. We also note that complicated 
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flavour structures that only appear in the — I h coefficients vanish in the J\f = 2 theory, 



n/(3 - Usf - 2(4 - n/)2 ^ (A/" - 2){M - 4f 
2(1 - 2ns) + nf -3{Af - 2). 



(3.15) 
(3.16) 



This shows that the integral basis for A/" = 2 is simpler than that of Af = 1 , a feature which 
is new for four-dimensional maximal cuts at three-loops. At one-loop both M = 1 and 
J\f = 2 have the same integral basis with no rational terms. At two- loops the two theories 
could differ in the lower propagator integrals but not in the maximal cut terms. We have 
also checked that by taking two-particle cuts the full integrand factorizes onto the two-loop 
results for the full integrand [11]. 

4 Conclusions 

In this paper we have shown that an integrand reduction technique based on computational 
algebraic geometry techniques can be generalized to three-loop amplitudes. As a first step 
in this direction we considered maximal cuts in four dimensions for the planar triple box 
topology contributing to massless 2 — > 2 scattering. 

By using primary decomposition the 14 independent branches of the 10 on-shell con- 
ditions could be found automatically. An explicit parametrisation of these solutions was 
found such that the integrand would take the form of a simple Laurent expansion in two free 
variables. The integrand on each solution factorises into a product of ten tree amplitudes 
which can be used to extract a set of 622 terms in the Laurent expansion. 

Having derived a minimal parametrisation of the integrand in terms of 7 irreducible 
scalar products we were able to construct an invertible matrix to map the coefficients of 
the Laurent expansion to the coefficients of the ISPs. Though the matrix was quite large 
it was possible to invert the system efficiently by using a branch-by-branch reconstruction 
using Grobner bases and the intersection of ideals. 

Finally we were able to reduce the whole expression to master integrals using the 
Reduze2 package for integration by parts identities. This whole procedure derives a com- 
plete reduction for an arbitrary process valid in any renormalizable gauge theory. As an 
application of the technique we computed the MI coefficients for gluon-gluon scattering in 
Yang-Mills theory with adjoint scalars and fermions. 

Though a long way from a complete four-point amplitude in QCD the computation 
presented here can be seen as a small step in the right direction. We hope the techniques 
presented here will be useful when making the necessary generalizations to L>-dimensional 
cuts and fewer particle cuts. 
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A Branch-by-branch polynomial fitting method 



Consider generalized unitarity for an L-loop diagram. Let R be the polynomial ring of 
ISPs and I be the ideal generated by cut equations. 

By the integrand-level reduction algorithm, the original integrand numerator is a poly- 
nomial N € R, while [N] is its image in the quotient ring R/I. [N] can be expanded over 
the integrand basis. The reduced integrand A^, which is the simplest representative for [N], 
{[N] = [N] in R/I), can be obtained by dividing A'' towards the Grobner basis G{I) of / 
[15, 16] . 

Often the cut solution have several branches. In other words, by primary decomposi- 
tion, I is decomposed to the intersection of several primary ideals [15], 

n 

I = f|/, (A.l) 

1=1 

where n is the number of branches. 

Note that / C /j, so /j contains more constraints than I itself. Hence the integrand 
can be reduced further on each branch. We have the map, 

p:R/I^ R/h © ... © R/In, 

[N] ^ ([iV]i,...,[iV]„), (A.2) 

where [. . .] and [. . .]i stand for the equivalence relations in R/I and R/Ii. If A^ is known, it 
is straightforward to get the simplest representative Ni for [N]i, by calculating the Grobner 
basis for each /j. In general, A'j contains much fewer terms than A^. 

So it is much easier to fit each Ni than A^, from the product of tree amplitudes. After 
all Ni are fitted, the goal is to determine A^ from A'^j's. The existence of A^ is guaranteed by 
the existence of the original numerator A^. We just need the following uniqueness condition: 

Lemma 1. The map p is injective. In other words, if N exists, then Ni, . . . ,Nn uniquely 
determine [N]. 

Proof. Assume that there are two polynomials A^ and A^' such that p{N) = p{N'). So 
[N]i = [N']i = [Ni]i, Mi. This means that A^ - A^' G U, Mi and 

n 

N -N' e{^Ii = I, (A.3) 

1=1 

so [N] = [N'] in R/I. □ 

This is an analogy of Chinese remainder theorem, however we do not need the coprime 
condition since the existence of A^ is guaranteed by physics. After [A^] is determined, it is 
straightforward to find its simplest representative A^ by Grobner basis G{I) of /. 

In practice, we present the following algorithm to get A^ from A^^i, . . . ,Nn. First, we 
consider the case n = 2, 
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1. For two ideals Ii = (/i, . . . , /m^) and I2 = {hi . . . hm2), Calculate the Grobner basis 
G{Ii + I2) = {gi, ■ ■ ■ gm} for the ideal Ii + l2- Record the transform matrix, 

mi m2 

9i = ^ o-ijfj + bikhk- (A.4) 
j=i k=i 

2. Divide A^i — N2 towards G{Ii + I2), where r is the remainder, 

Nl-N2 = J2^ig^ + r. (A.5) 

3. If r = 0, rewrite 

i j i k 

where F\ G /i and F2 (z l2- Then N = Ni — Fi. If r 7^ 0, print warning message and 
stop. 

4. Calculate the intersection Ii n/2 and its Grobner basis G{Ii n/2). Divide N towards 
G{Ii n I2) and the remainder is A^. 

The validity of this algorithm is self-evident since N = Ni — Fi = N2 + F2 so [N]i = [Ni]i 
and [N]2 = [A^2]2- And as long as A^ exists, r must be zero. 

For cases with more than 2 branches, we just need to repeat the above algorithm for 
n — 1 times, 

1. Let J = Ii, / = A^i. 

2. For z = 1 to n — 1 

• Carry out the 2-branch algorithm for polynomials (/, Aj+i) and the two ideals 
(J, Ij+i). Then redefine / as the output polynomial. 

• J = Jn/j+i. 

3. A = /. 

The validity can be checked by induction. We realise this algorithm in a Macaulay2 program 
[28]. 

B Flavour Configurations in the Triple Box 

The pre-factors of the 34 configurations appearing in (3.1) are given in Table 4. We note 
that by combining these configurations together with different colour factors the results 
presented here would also be valid in QCD. The symmetry factor of —1 for each fermion 
loop is included in the pre- factor. 
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Table 4. The definitions of the flavour pre- factors in eq. (3.1). 



For completeness we also present the well-known tree level amplitudes used in this 
paper, 

(12)3 



— ijdg (Ig, 2s, 3g ) 

l(0) 



(23) (31) 

(13)! 
(12) 

(13) (23) 
(12) 



i4"^(l-,2+,3,) = (12) 



(B.l) 

(B.2) 

(B.3) 
(B.4) 



-i^4 \lg , 2^, 3^, 4g ) 
-i^4 \lg , 2^, 3g , 4^) 



(12)= 



(23) (34) (41) 

(41)^ 
(12)(23)(34) 

(13)^ 
(12)(23)(34)(41) 



(B.5) 
(B.6) 
(B.7) 



References 

[1] Z. Bern, J. Carrasco, L. J. Dixon, H. Johansson, and R. Roiban, The Ultraviolet Behavior of 
N=8 Supergravity at Four Loops, Phys. Rev. Lett. 103 (2009) 081301, [arXiv: 0905. 2326]. 



- 12 - 



[2] Z. Bern, J. Carrasco, L. J. Dixon, H. Johansson, and R. Roiban. The Complete Four-Loop 
Four-Point Amplitude in N=4 Super-Yang-Mills Theory, Phys.Rev. D82 (2010) 125040, 
[arXiv: 1008. 3327]. 

[3] Z. Bern, S. Davies, T. Dennen, and Y.-t. Huang, Absence of Three-Loop Four-Point 

Divergences in N^4 Supergravity, Phys.Rev.Lett. 108 (2012) 201301, [arXiv: 1202.3423]. 

[4] J. J. M. Carrasco and H. Johansson, Generic multiloop methods and application to N=4 
super- Yang- Mills, J.Phys.A A44 (2011) 454004, [arXiv : 1103 . 3298]. 

[5] Z. Bern, L. J. Dixon, D. C. Dunbar, and D. A. Kosower, One loop n point gauge theory 
amplitudes, unitarity and collinear limits, Nucl.Phys. B425 (1994) 217-260, 
[hep-ph/9403226]. 

[6] Z. Bern, L. J. Dixon, D. C. Dunbar, and D. A. Kosower, Fusing gauge theory tree amplitudes 
into loop amplitudes, Nucl.Phys. B435 (1995) 59-101, [hep-ph/9409265]. 

[7] R. Britto, F. Cachazo, and B. Feng, Generalized unitarity and one-loop amplitudes in N=4 
super- Yang- Mills, Nucl.Phys. B725 (2005) 275-305, [hep-th/0412103]. 

[8] G. Ossola, C. G. Papadopoulos, and R. Pittau, Reducing full one-loop amplitudes to scalar 
integrals at the integrand level, Nucl.Phys. B763 (2007) 147-169, [hep-ph/0609007]. 

[9] P. Mastroha and G. Ossola, On the Integrand- Reduction Method for Two-Loop Scattering 
Amplitudes, JHEP 1111 (2011) 014, [arXiv : 1107 . 6041]. 

[10] D. A. Kosower and K. J. Larsen, Maximal Unitarity at Two Loops, arXiv: 1108. 1180. 

[11] S. Badger, H. Frehesvig, and Y. Zhang, Hepta-Guts of Two-Loop Scattering Amplitudes, 
JHEP 1204 (2012) 055, [arXiv : 1202 . 2019]. 

[12] K. J. Larsen, Global Poles of the Two-Loop Six-Point N=4 SYM integrand, 
arXiv: 1205.0297. 

[13] S. Caron-Huot and K. J. Larsen, Uniqueness of two -loop master contours, arXiv : 1205 . 0801. 

[14] R. H. Kleiss, I. Malamos, C. G. Papadopoulos, and R. Verheyen, Counting to one: 
reducihility of one- and two-loop amplitudes at the integrand level, arXiv: 1206.4180. 

[15] Y. Zhang, Integrand- Level Reduction of Loop Amplitudes by Computational Algebraic 
Geometry Methods, arXiv : 1205 . 5707. 

[16] P. Mastrolia, E. Mirabella, G. Ossola, and T. Peraro, Scattering Amplitudes from 
Multivariate Polynomial Division, arXiv : 1205 . 7087. 

[17] K. Chetyrkin and F. Tkachov, Integration by Parts: The Algorithm to Calculate beta 
Functions in 4 Loops, Nucl.Phys. B192 (1981) 159-204. 

[18] T. Gehrmann and E. Reniiddi, Differential equations for two loop four point functions, 
Nucl.Phys. B580 (2000) 485-518, [hep-ph/9912329]. 

[19] J. Gluza, K. Kajda, and D. A. Kosower, Towards a Basis for Planar Two-Loop Integrals, 
Phys.Rev. D83 (2011) 045012, [arXiv : 1009 . 0472]. 

[20] R. M. Schabinger, A New Algorithm For The Generation Of Unitarity- Compatible 
Integration By Parts Relations, JHEP 1201 (2012) 077, [arXiv : 1111 . 4220]. 

[21] S. Laporta, High precision calculation of multiloop Feynman integrals by difference equations, 
Int.J.Mod.Phys. A15 (2000) 5087-5159, [hep-ph/0102033]. 



-13- 



[22] C. Anastasiou and A. Lazopoulos, Automatic integral reduction for higher order perturbative 
calculations, JHEP 0407 (2004) 046, [hep-ph/0404258]. 

[23] A. Smirnov, Algorithm FIRE - Feynman Integral REduction, JHEP 0810 (2008) 107, 
[arXiv: 0807. 3243]. 

[24] C. Studerus, Reduze- Feynman Integral Reduction in C++, Comput.Phys.Commun. 181 
(2010) 1293-1300, [arXiv:0912.2546]. 

[25] A. von Manteuffel and C. Studerus, Reduze 2 - Distributed Feynman Integral Reduction, 
arXiv: 1201.4330. 

[26] C. Bcrgcr, Z. Bern, L. Dixon, F. Febres Cordero, D. Forde, et. ai. An Automated 

Implementation of On-Shell Methods for One-Loop Amplitudes, Phys.Rev. D78 (2008) 
036003, [arXiv: 0803. 4180]. 

[27] P. Mastrolia, G. Ossola, C. Papadopoulos, and R. Pittau, Optimizing the Reduction of 
One-Loop Amplitudes, JHEP 0806 (2008) 030, [arXiv : 0803 . 3964]. 

[28] D. R. Grayson and M. E. Stillman, "Macaulay2, a software system for research in algebraic 
geometry." Available at http://www.math.uiuc.edu/Macaulay2/. 

[29] V. A. Smirnov, Analytical result for dimensionally regularized massless on shell planar triple 
box, Phys.Lett. B567 (2003) 193-199, [hep-ph/0305142]. 

[30] Z. Bern, L. J. Dixon, and V. A. Smirnov, Iteration of planar amplitudes in maximally 

supersymmetric Yang-Mills theory at three loops and beyond, Phys.Rev. D72 (2005) 085001, 
[hep-th/0505205]. 



-14- 



